GASTO ÓPTIMO EN RIEGO POR MELGAS EN 
PRESENCIA DE UN MANTO FREÁTICO SOMERO 


» Heber Saucedo* e 
Instituto Mexicano de Tecnología del Agua 


*Autor de correspondencia 


e Manuel Zavala e 
Universidad Autónoma de Zacatecas, México 


e Carlos Fuentes + 
Universidad Autónoma de Querétaro, México 


Resumen 


El gasto óptimo de riego, es decir, aquel para el cual se obtiene la mayor eficiencia 
de uniformidad manteniendo valores elevados de las eficiencias de aplicación y 
de requerimiento de riego, se determina para el caso en que se tiene la presencia 
de un manto freático en el perfil de suelo en riego por melgas. El cálculo se efectúa 
haciendo uso de un modelo de simulación desarrollado por los autores de este 
trabajo, el cual hace posible describir las cuatro fases del riego por melgas, y 
tiene como base el acoplamiento numérico de las ecuaciones de Saint-Venant y la 
ecuación de Richards; las primeras describen el flujo del agua con una superficie 
libre que tiene lugar sobre el suelo y la segunda permite simular la transferencia de 
agua en el suelo. Se ha obtenido que no obstante la presencia de un manto freático 
en el perfil del suelo, la relación entre la longitud de la melga y el gasto óptimo es 
básicamente lineal, situación que se presenta en el caso de un perfil de suelo semi- 
infinito; en otras palabras, se deduce que la dependencia básicamente lineal que 
guardan la longitud de la melga y el gasto óptimo de riego es independiente de la 
condición inicial de distribución de las presiones en el suelo. 


Palabras clave: riego por melgas, ecuaciones de Saint-Venant, ecuación de Richards. 


Introducción 


La descripción del flujo del agua en el riego por 
melgas requiere tener en cuenta la ecuación de 
Richards para describir el flujo del agua en el 
suelo y las ecuaciones de Saint-Venant para el 
flujo del agua sobre la superficie del suelo, de 
tal forma que su acoplamiento permita conocer 
las formas de los perfiles de flujo superficial y 
la distribución de los potenciales de presión 
en el suelo según avanza el riego. El tirante 
proporcionado por las ecuaciones de Saint- 
Venant se utiliza para definir una condición 
de frontera tipo Dirichlet para la ecuación 
de Richards unidimensional, cuya solución 
numérica permite el cálculo de la lámina 
infiltrada necesaria para la solución de las 


ecuaciones de Saint-Venant, teniéndose por lo 
tanto un procedimiento iterativo. 

Saucedo et al. (2003) han analizado el efecto 
que tiene la posición inicial de un manto freático 
sobre la evolución del frente de avance en el 
riego por melgas, concluyendo que dicho efecto 
es mayor cuanto más cercana sea la posición 
inicial del manto freático en relación con la 
superficie del suelo. Una extensión natural de 
dicho análisis es calcular el efecto que tiene 
la posición inicial del manto freático no sólo 
en el frente de avance, sino en las cuatro fases 
que constituyen el riego por melgas y, en 
consecuencia, en la eficiencia del riego. 

El objetivo del presente trabajo es utilizar 
el acoplamiento de las ecuaciones unidimen- 
sionales de Saint-Venant y Richards para calcular 
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el gasto óptimo de riego que debe aplicarse 
cuando se tiene la presencia de un manto 
freático somero en el perfil del suelo, es decir, 
según la posición inicial del manto freático y la 
lámina de riego por aplicar. 


Flujo del agua sobre la superficie del 
suelo 


El flujo del agua con una superficie libre 
es modelado con las ecuaciones de Saint- 
Venant, que resultan de la aplicación de las 
leyes de conservación de masa y cantidad 
de movimiento. En una melga, la relación 
entre su ancho y el tirante de agua permite 
considerar las ecuaciones correspondientes al 
escurrimiento sobre una superficie de ancho 
infinito (Woolhiser, 1975). La ecuación de 
continuidad se escribe como: 


0 (1) 
ot 9x ot 


La ecuación de momentum se escribe en la 
forma recomendada por Saucedo et al. (2005): 


log 2q0 2 ah 
1 


hot ox 1 

q ol 
+2 =0 2 
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donde ql(x,t) = Ul(x,t)h(x,t) es el gasto por 
unidad de ancho de melga [L?T”"]; x, la 
coordenada espacial en la dirección principal 
del movimiento del agua en la melga [L]; t, el 
tiempo [T]; U, la velocidad media; h, el tirante 
de agua [L]; J,, la pendiente topográfica de la 
melga [LL”]; J, la pendiente de fricción [LL”); 
V, =01/0t, el flujo de infiltración [LT”], es decir, 
el volumen de agua infiltrado en la unidad de 
tiempo por unidad de ancho y por unidad 
de longitud de la melga; 1, la lámina infiltra- 
da [Ll; g, la aceleración gravitacional [LT]; 
p = 1 - a, el parámetro adimensional, siendo 
a = 1 - U,,/U, donde U,, es la proyección en 


la dirección del movimiento de la velocidad 
de salida de la masa de agua debido a la 
infiltración. 

La relación entre las variables hidráulicas q 
y h con la pendiente de fricción, denominada 
ley de resistencia hidráulica, es adoptada de 
acuerdo con Fuentes et al. (2004), haciendo uso 
de una ley potencial de resistencia: 


3 d 
q> 18 6) 


donde v es el coeficiente de viscosidad 
cinemática del agua [L?T"] y k es un factor 
adimensional; d es un parámetro adimensional, 
que varía entre 1/2 < d < 1 en función del tipo 
de flujo; los valores extremos d = 1/2 y d = 1 
corresponden, respectivamente, al régimen de 
Chézy y al flujo laminar de Poiseuille. 

Para establecer las condiciones inicial y de 
frontera que deben sujetar a las ecuaciones de 
Saint-Venant para modelar el riego por melgas, 
se tiene en cuenta una melga cerrada, para 
evitar la pérdida de masa fuera del dominio de 
riego. Para la fase de avance se tiene: 


g(x,0)=0 y  h(x,0)=0 (4) 


q(0,t)=3,, q(x,,t)=0 y h(xj,t)=0 (5) 


donde xt) es la posición del frente de onda 

para el tiempo t y q, es el gasto de aporte 

en la entrada de la melga. Para la fase de 
almacenamiento: 

q(0,t)=4,,  q(L,)=0 (6) 

donde L es la longitud de la melga y q, es el 

gasto de aporte en la cabecera de la melga. Para 

la fase de consumo: 

q(0,t)=0,  q(L,t)=0 (7) 

Para cerrar el sistema es necesario conocer 

la forma en que evoluciona en el tiempo la 
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lámina infiltrada en toda posición sobre la 
melga, es decir, la ley de infiltración, la cual 
es proporcionada con base en la solución de la 
ecuación de Richards (1931). 


Flujo del agua en el suelo 


La descripción del proceso de infiltración del 
agua en el suelo puede efectuarse teniendo 
como base la ecuación de Richards (1931), 
que resulta de la combinación de la ecuación 
de continuidad con el campo de velocidades 
calculado conforme a la ley de Darcy, y que en 
su forma tridimensional, sin tener en cuenta la 
extracción de agua por las plantas, se escribe 
como: 


oy dK dy 
C(w—É=VW-[K(w)Vy]-2 
(y) 0 [ (y) y] 


(8) 
dy 0z 


donde y es el potencial de presión del agua 
en el suelo, expresado como la altura de una 
columna equivalente de agua [L] (positivo en la 
zona saturada y negativo en la zona no saturada 
del suelo); C(y) = d0/dy es denominada la 
capacidad específica de humedad del suelo; 8 
= O(y) es el volumen de agua por unidad de 
volumen de suelo o contenido volumétrico de 
agua [L*L*], y es una función de y conocida 
como curva característica de humedad o 
curva de retención del agua; K = K(y) es la 
conductividad hidráulica [LT“], que en un 
suelo parcialmente saturado es una función del 
potencial de presión; el potencial gravitacional 
es asimilado a la coordenada espacial z, 
orientada positivamente hacia abajo [L]; V = 
(9/9x,0/0y,0/0z) es el operador gradiente; x, y 
son las otras dos coordenadas espaciales [L] ,y 
tes el tiempo [T]. 

El riego es un fenómeno que se efectúa en 
tres dimensiones, lo cual hace necesario que el 
flujo del agua en el suelo sea descrito utilizando 
la ecuación (8), sin embargo, es conveniente 
aceptar la hipótesis de que el fenómeno se 


efectúa en planos paralelos al desarrollo de la 
melga, y entonces es posible utilizar la forma 
bidimensional de la ecuación de Richards: 


y 0 ayl a oy 
Cde? Ki z | 2 0 o i] 0) 
Habida cuenta de los resultados de Saucedo 
et al. (2006) sobre el grado de aproximación 
de la hipótesis del tiempo de contacto en el 
riego por melgas, es posible simplificar más 
la aproximación, empleando la ecuación de 
Richards en su forma unidimensional: 


dy_0 dp _ 
AE 1] (10) 


que debe resolverse sobre el dominio de 
solución definido por una columna de suelo. 

Como condición inicial para la solución de 
la ecuación de Richards unidimensional, se 
debe especificar la distribución de las presiones 
en el espacio; se considera una distribución 
hidrostática de las presiones, teniendo como 
referencia la posición del manto freático, es 
decir: 

y=z-P, Osz<Ec, t>0 (11) 
donde P es la profundidad del manto freático y 
Ec es el espesor de la columna de suelo. 

En la superficie del suelo se impone una 
condición de frontera tipo Dirichlet, con el 
potencial de presión igual al tirante de agua 
calculado mediante las ecuaciones de Saint- 
Venant: 
t>0 (12) 

La condición de frontera inferior se asigna 
como condición de Neumann con flujo nulo: 


t>0 


z=Ec, (13) 
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teniendo en cuenta que el espesor de la columna 
(Ec) debe ser mayor que la profundidad del 
manto freático. 

La solución de la ecuación de Richards hace 
indispensable representar las propiedades 
hidrodinámicas del suelo, expresando el 
potencial de presión (yw) como una función 
del contenido volumétrico de agua (0) y la 
conductividad hidráulica K como una función 
de 6. Como es señalado por Fuentes et al. 
(1992), la combinación de las características 
hidrodinámicas de Fujita (1952) y Parlange et 
al. (1982), es conveniente en estudios teóricos, 
como la construcción de soluciones analíticas 
exactas, y en estudios experimentales, el uso 
de la combinación de la curva de retención 
propuesta por (1980), 
considerando la restricción de Burdine (1953) 
con la curva de conductividad hidráulica 
propuesta por Brooks y Corey (1964), debido a 


van  Genuchten 


que satisfacen las propiedades integrales de la 
infiltración y a la facilidad para la identificación 
de sus parámetros. 

La curva de retención propuesta por 
van Genuchten (1980) 


[o(1)-0,]/[9,-0,]= 


es un valor característico de la presión del agua 


se escribe como 


1 +(1y/ w)']” donde yv, 


en el suelo; 1 y n son dos parámetros de forma 
empíricos, relacionados por la restricción de 
Burdine (1953): m = 1 -— 2/n, con0<m<l, y 
n > 2, 6, es el contenido volumétrico de agua 
a saturación efectiva del suelo y 0, es el conte- 
nido volumétrico de agua residual. La 
conductividad hidráulica propuesta por 
Brooks y Corey (1964) se representa como 
K(0)=K.[(9-0,)/(0, 0,31", donde n es un 
parámetro de forma empírico y positivo. 

La lámina infiltrada se calcula para cada 
posición x, donde se resuelven de forma 


numérica las ecuaciones de Saint-Venant: 
FE 

I(x,,4)= S[9(x,,2,t) -0(x,,2,0)]az, donde F es 
0 

la posición del frente de humedecimiento y 


9(x, z, 0) es el contenido volumétrico de agua 
inicial. 


Solución numérica 


Solución numérica de las ecuaciones de Saint- 
Venant 


Para el cálculo del gasto óptimo se usa un 
esquema numérico para las cuatro fases del 
riego por melgas desarrollado por Saucedo 
(2003). La forma discreta de la ecuación de 
continuidad para la fase de avance se escribe 
como: 


(09, +(1 -00)g, J9t-(x, -x,)[0(», +1,) 


+ (1-0) Mi +15 )] (14) 
- [0h, + 1-4), +91,+(1-4)1,](x, -x,) 
[9h +(1-)h,, + 91, +(1-0)1,n ](%n =x;) =0 


para las fases de almacenamiento, consumo y 
recesión se escribe como: 


(09, +(1-00)g,]ot -[09, +(1-00)g,, J9t 


- [9h, +(1-4)h, + bl, +(1-4)1, ](x, xy) 
+ | 9h, +(1-9)h,, + 91, +(1-9)1,, (xy =x;)=0 


(15) 


La ecuación de momentum guarda la misma 
forma discreta para las cuatro fases del riego: 


23h (9, -q,)0++ (gh 7), -h,)0t 
He (x, -x,)[o9, +(1-0)4, 09, -(1-0)4,] 
+90) (x, -x,)[0J, +(1-0)J,-J,] 


+B3+7h (x,-x,)[01,+(1-0)1,-01,,=(1-0)1,|=0 


(16) 


El cálculo de los coeficientes se realiza con 
base en los valores pertenecientes al nivel de 
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tiempo anterior: 7 = (1 — p),, + d9, h= (1 - p) 
h, + dh, (Saucedo et al., 2005). Se introducen 
las variables pequeñas (Strelkoff y Katopodes, 
1977): 6h, 9q,, 8h, y 0q,, tales que, para las celdas 
interiores, se tiene: h, = h, + 0h, h,=h,, + 0h, q,= 
q,+ 04, y q,=4,, + 0q,. Se sustituyen en las formas 


m 


discretas de las ecuaciones de continuidad 
y momentum, y se obtiene un sistema de 
ecuaciones algebraicas, cuya solución permite 
avanzar en el tiempo los valores de las variables 
hidráulicas. 


Solución numérica de la ecuación de Richards 


La ecuación de Richards unidimensional (10) 
es discretizada en el espacio utilizando ele- 
mentos finitos y en el tiempo mediante un 
esquema implícito en diferencias finitas. El 
procedimiento se encuentra bien documentado 
en la literatura y puede consultarse, por 
ejemplo, en Saucedo et al. (2002). 


Aplicaciones 


Caracterización hidrodinámica y 
discretización 


El efecto de la presencia de un manto freático 
en el cálculo del gasto óptimo en el riego por 
melgas se ilustra para el caso de la zona de 
riego de La Chontalpa, Tabasco, México; se 
hace uso de la información reportada para 
un experimento realizado por Pacheco (1995): 
gasto total de aporte en la entrada Q,=15 l/s; 
ancho de la melga A = 10.5 m; gasto unitario 
de aporte en la entrada q, = 0.001428 m*/s/m; 
longitud de la melga L = 100 m; pendiente de 
la melga J, = 0.00085 cm/cm. 

La caracterización hidrodinámica de la 
arcilla de La Chontalpa ha sido efectuada por 
Saucedo et al. (2003). Se tienen los valores 
siguientes para la característica de humedad de 
van Genuchten (1980): contenido volumétrico 
de agua residual 6, = 0.0 cm? /cm?; contenido 
volumétrico de agua a saturación 0, = 0.5245 
0.066. 
Para la ecuación de Brooks y Corey se tiene 


cm?/cm?; parámetro de forma m = 


el parámetro n = 21.19. Los dos parámetros 
restantes (y, K.) se estiman a fin de reproducir 
la curva de avance observada de modo 
experimental en una prueba de riego. Se ha 
utilizado un valor f =2 del parámetro asociado 
con la extracción de momentum debida a la 
infiltración. El cálculo del número de Reynolds 
permite clasificar el régimen de flujo laminar, 
por lo cual se utiliza d = 1 como parámetro de 
la ley de resistencia. 

La suma de los cuadrados de los errores 
entre la curva de avance calculada con el 
sistema acoplado de Saint-Venant y Richards, 
y la curva de avance observada es mínima con 
los valores y», =-15 cm y K_ = 1.86 cm/h. En la 
figura 1 se observa la buena concordancia entre 
las dos curvas de avance, R?= 0.9984, obtenida 
utilizando un coeficiente adimensional de 
la ley de resistencia k = 1/300. El detalle del 
procedimiento se puede consultar en Saucedo 
et al. (2003). 

Se ha utilizado un paso de tiempo constante 
0t = 1.0 s tanto para las ecuaciones de Saint- 
Venant como para la ecuación de Richards. 
Para la ecuación de Richards se utilizaron pasos 
. = 0.02 cm 


mín 


= 1.0 cm. La discretización utilizada 


de espacio mínimo y máximo: Az 
y Az 
guarda semejanza con las reportadas en la 
literatura (e.g., Katopodes y Strelkoff, 1977): 
Ot... = 5 s; Akanbi y Katopodes (1988): 0t_. 
= 1 s, para las cuales se reportan resultados 


máx 


XxX 


consistentes en las variables hidráulicas del 
riego. La discretización para la ecuación de 
Richards es más fina que la reportada por 
Saucedo et al. (2002), quienes comparan perfiles 
de humedad proporcionados por una solución 
analítica exacta con los obtenidos mediante 
una solución numérica para la ecuación de 
Richards, los valores del error relativo máximo 
para el caso de dos suelos contrastantes; una 
arcilla clara de Yolo, Estados Unidos, y una 
arena del río Isére, Francia, presentan una cota 
superior de 0.25% cuando se utiliza un paso de 
tiempo máximo 0f, . =30 s y pasos de espacio 
mínimo y máximo Az, = 0.04 cm y Az, . = 
10.0 cm. El suelo arcilla de la Choltalpa tiene 
un comportamiento hidrodinámico similar al 
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Figura 1. Calibración del modelo; y, =-15 cm y k, =1.86 cm/h, R?=0.9984. 


de la arcilla clara de Yolo; en consecuencia, con 
la discretización utilizada, el error relativo 
máximo debe tener como una cota superior 
0.25%. 


Efecto de la presencia de un manto freático en 
riego por melgas 


La representación de la presencia de un manto 
freático somero con la condición de frontera 
dada por la ecuación (13), hace que en todo 
punto ubicado sobre la melga solamente sea 
posible que se tenga una lámina infiltrada 
máxima, definida en relación con la diferencia 
que existe entre la distribución inicial del 
contenido volumétrico de agua en el suelo y 
un perfil de suelo saturado. Cuando se alcanza 
esta lámina máxima se presenta un aumento 
en la velocidad de la onda del flujo a superficie 
libre; según puede apreciarse en la figura 2, 
se ha tenido en cuenta el suelo arcilla de La 
Chontalpa, con una profundidad del manto 
freático de 100 cm. 


Relación entre la longitud de la melga y el 
gasto óptimo 


El acoplamiento de las ecuaciones de Saint- 
Venant y Richards en el riego por melgas es 
utilizado con la finalidad de calcular el gasto 
óptimo de riego, es decir, aquel para el que se 
obtiene el mayor coeficiente de uniformidad, 
manteniendo valores aceptables de las 
eficiencias de aplicación y de requerimiento de 
riego, cuando se tiene la presencia de un manto 
freático en el perfil del suelo. La eficiencia de 
uniformidad ha sido evaluada utilizando el 


coeficiente de uniformidad de Christiansen 
(CUC): CUC=1-(1/m1))/1,- 1 donde 1, es la 
¡=1 


lámina infiltrada en el punto ¡[L]; 7 es la lámina 
infiltrada media [L], y n es el número de puntos 
considerados para efectuar el cálculo. 

En la figura 3 se muestra que la relación 
entre el gasto óptimo y la longitud de la 
melga para la cual es calculado dicho gasto es 
suficientemente bien descrita utilizando una 
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Tirante (cm) 


Lámina infiltrada (cm) 


Distancia (cm) 


Figura 2. Efecto de la presencia de un manto freático somero en los perfiles de flujo de un riego por melgas. Perfiles de flujo 


cada diez minutos. Suelo arcilla de La Chontalpa; profundidad del manto freático, 100 cm. 


correspondencia lineal; la figura se ha hecho 
teniendo en cuenta una lámina de riego de 10 
cm, pero el comportamiento es similar, con in- 
dependencia de la lámina de riego que se desea 
aplicar. El resultado es similar al obtenido 
para el caso en que no se presenta un manto 
freático en el perfil del suelo, lo cual indica que 
la dependencia básicamente lineal que guarda 
la longitud de la melga y el gasto óptimo de 
riego es independiente de la condición inicial 
de distribución de las presiones en el suelo. 
A manera de ejemplo, puede apreciarse en 
la figura 4 la distribución final de la lámina 
infiltrada para una lámina de riego de 10 cm 
y una profundidad de manto freático de 100 
cm en el suelo arcilla de La Chontalpa, dicha 
profundidad ha sido usada para definir la 
distribución inicial de las presiones en el suelo 
conforme a la ecuación (11). 


Cuadro de diseño del riego por melgas 
para un suelo arcilloso con presencia de 
manto freático 


En el cuadro 1 se muestran los valores del gasto 
unitario; es decir, gasto que debe aplicarse por 
unidad de ancho y por unidad de longitud 
de la melga, correspondientes a diferentes 
láminas de riego y profundidades iniciales del 
manto freático, las casillas vacías indican que la 
cantidad de agua que se desea aplicar es mayor 
que la cantidad de agua que puede infiltrarse en 
el suelo. En el cuadro 2 se muestran los tiempos 
de riego correspondientes a los gastos unitarios 
mostrados en el cuadro 1. Las posiciones 
iniciales del nivel freático consideradas son 
100, 150, 200 y 10 000 cm, y puede apreciarse 
que existe una tendencia creciente del gasto 
óptimo conforme se incrementa la profundidad 
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Figura 3. Relación entre la longitud de melga y el gasto de aporte óptimo para el suelo arcilla de La Chontalpa para 
una lámina de riego de 10 cm y una profundidad inicial del manto freático de 200 cm; K_ se expresa en m/s. 
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Figura 4. Distribución final de la lámina infiltrada en el suelo arcilla de La Chontalpa. Profundidad del manto freático, 
100 cm; lámina de riego, 10 cm; CUC = 0.96. 
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Cuadro 1. Gasto de riego (1/s/m) para la aplicación óptima de la lámina de riego. 


PNF (cm) Lr =3 Lr=5 Lr=6 Lr =8 Lr =10 Lr=12 Lr =14 
100 0.01162 0.01083 ii * > 7 ij 
150 0.01455 0.01110 0.00896 0.00607 0.00441 ES ES 
200 0.01712 0.01141 0.00927 0.00688 0.00524 0.00445 0.00388 
10 000 0.03236 0.01976 0.01554 0.01221 0.01001 0.00826 0.00775 
Cuadro 2. Tiempo de riego (h) para la aplicación óptima de la lámina de riego. 
PNF (cm) Lr =3 Lr=5 Lr=6 Lr =8 Lr=10 Lr=12 Lr=14 
100 0.450 1.445 * id * * id 
150 0.312 1.389 2.067 4.073 6.862 E S 
200 0.289 1.339 2.050 3.478 5.884 8.112 11.039 
10 000 0.022 0.817 1.200 DAME) 3.012 4.578 5.505 


del manto freático. Debido a que las pro- 
fundidades consideradas no se encuentran en 
la zona asintótica de la curva de retención de 
humedad, que se forma desde el manto freático 
hacia la superficie del suelo como consecuencia 
de la distribución de presiones hidrostática 
adoptada, no es posible alcanzar el previsible 
comportamiento asintótico del gasto óptimo 
de riego. No obstante lo anterior, los valores 
obtenidos pueden servir como guía para el 
manejo del agua en la zona productora de 
arroz de La Chontalpa, teniendo en cuenta las 
hipótesis básicas de que el suelo es homogéneo 
y que constituye un medio indeformable. 


Conclusiones 


El gasto óptimo de riego, es decir, aquel 
para el cual se obtiene la mayor eficiencia 
de uniformidad, manteniendo valores acep- 
tables de las eficiencias de aplicación y de 
requerimiento de riego, se determina para el 
caso en que se tiene la presencia de un manto 
freático en el perfil de suelo en riego por 
melgas. El cálculo se efectúa haciendo uso de 
un modelo de simulación desarrollado por los 
autores, el cual hace posible describir las cuatro 
fases del riego por melgas, y tiene como base 
el acoplamiento numérico de las ecuaciones 


de Saint-Venant y la ecuación de Richards; las 
primeras describen el flujo del agua con una 
superficie libre que tiene lugar sobre el suelo 
y la segunda permite simular la transferencia 
de agua en el suelo. Se ha obtenido que, no 
obstante la presencia de un manto freático en 
el perfil del suelo, la relación entre la longitud 
de la melga y el gasto óptimo es básicamente 
lineal, situación que se presenta en el caso de un 
perfil de suelo semi-infinito; en otras palabras, 
se deduce que la dependencia básicamente 
lineal que guardan la longitud de la melga y 
el gasto óptimo de riego es independiente 
de la condición inicial de distribución de las 
presiones en el suelo. 
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Abstract 


SAUCEDO, H., ZAVALA, M. € FUENTES, C. Optimal flow for border irrigation with the 
presence ofa shallow water table. Water Technology and Sciences (in Spanish). Vol. IV, No. 
5, November-December, 2013, pp. 99-109. 


Optimal irrigation flow-—1.e. the flow resulting in the greatest uniformity efficiency while 
maintaining acceptable application efficiency and irrigation requirement values-was 
determined for border irrigation with a water table in the soil profile. The calculation was 
performed using a simulation model developed by the authors of this work to describe the 
four border irrigation phases, and was based on the coupling of Saint-Venant and Richards 
equations. The former equation describes the flow of water with a free surface over the soil and 
the latter enables simulating the transfer of water into the soil. Regardless of the presence of a 
water table in the soil profile, the relationship between the border length and the optimal flow 
is basically linear, a situation which occurs when the soil profile is semi-infinite. That is, it is 
deduced that the basically linear dependence between the border length and optimal irrigation 
flow is independent of the initial distribution of the pressure in the soil. 


Keywords: border irrigation, optimal irrigation flow, Richards equation, Saint-Venant 
equations, shallow water table. 
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